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^-H We investigate standard and non-standard solar neutrino signals in direct dark matter 

detection experiments. It is well known that even without new physics, scattering of solar 
neutrinos on nuclei or electrons is an irreducible background for direct dark matter searches, 

O once these experiments reach the ton scale. Here, we entertain the possibility that neutrino 

interactions are enhanced by new physics, such as new light force carriers (for instance a "dark 
photon") or neutrino magnetic moments. We consider models with only the three standard 
neutrino flavors, as well as scenarios with extra sterile neutrinos. We find that low-energy 
neutrino-electron and neutrino-nucleus scattering rates can be enhanced by several orders 
of magnitude, potentially enough to explain the event excesses observed in CoGeNT and 

f~| CRESST. We also investigate temporal modulation in these neutrino signals, which can arise 

from geometric effects, oscillation physics, non-standard neutrino energy loss, and direction- 

^ ^ dependent detection efficiencies. We emphasize that, in addition to providing potential 

explanations for existing signals, models featuring new physics in the neutrino sector can 
i-S^ also be very relevant to future dark matter searches, where, on the one hand, they can be 

probed and constrained, but on the other hand, their signatures could also be confused with 
y—i dark matter signals. 

> 

PACS numbers: 

o 

\q 

^ 1. INTRODUCTION AND MOTIVATION 

(N 

Experiments searching for dark matter and those studying solar neutrinos have coexisted for 
many years, and even though many of the technological challenges — such as the suppression of 
^ radioactive backgrounds and the lowering of the energy threshold — are similar for both types of 

^ experiments, their physics programs have had little overlap. In solar neutrino physics, the size 

of the detector is more important than the energy threshold, whereas the direct search for dark 
matter scattering is only possible for low energy thresholds of order 10 keV or below. Until very 
recently, achieving such a low threshold has only been possible in small detectors with masses of 
only a few kg, too small to be of interest to neutrino physics. However, as dark matter detectors 
are becoming more and more massive, this is about to change. In ton scale experiments, solar 
neutrinos will constitute a non-negligible irreducible background p] which can only be overcome 
in detectors with directional sensitivity. 

In this work we consider the additional possibility that there is new physics in the neutrino 
sector, and we will examine how this may affect the signals observed by present and future dark 
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matter experiments. In particular, we will discuss models featuring a new light 1 GeV) gauge 
boson A' . We will also discuss the possible existence of neutrino magnetic moments and of extra 
("sterile") neutrino species coupled to the new gauge boson. ^ Sterile neutrinos are motivated by 
a number of anomalous results from short-baseline neutrino oscillation experiments and by 

cosmology |12H14j . whereas new light force carriers appear for instance in models of Sommerfeld- 
enhanced dark matter annihilation (see e.g. |15H19j ). Both sterile neutrinos and light force carriers 
were used by Nelson and Walsh to explain some of the neutrino oscillations anomalies |20l I21j 
(see also |22| for related work). As also pointed out by Pospelov [23], a new light force carrier 
can enhance the interaction rates of active and sterile neutrinos at low energy by several orders 
of magnitude compared to Standard Model neutrinos. Additionally, it can have a large impact on 
the Mikheyev-Smirnov-Wolfenstein (MSW) potential that neutrinos feel when propagating through 
matter, thus affecting oscillations. 

As examples, we will consider potential explanations of current dark matter anomalies, in 
particular the event excesses observed by CoGeNT [251 ES] and CRESST [26], as well as the 
annual modulation signals observed in CoGeNT [2S] and DAMA [27j. While we will see that 
finding a model that is consistent with all current direct detection data and with constraints from 
neutrino oscillation experiments is difficult (a fate that our scenario shares with most dark matter 
models), we will also see that future experiments can be considerably affected by new physics in 
the neutrino sector. It will be important for us that dedicated neutrino experiments have a much 
higher energy threshold than dark matter detectors, which can detect electron recoil energies as 
low as 0.5 keVee and nuclear recoil energies down to 2 keVnr^ in the case of CoGeNT [21] and 
CDMS [283. Thus, if new physics in the neutrino sector manifests itself most strongly at low energy, 
neutrino detectors would be insensitive, while dark matter detectors have excellent prospects of 
discovering or constraining such models. 

The outline of the paper is as follows: In section [2| we review neutrino scattering in the Stan- 
dard Model and compute the expected neutrino-electron and neutrino-nucleus scattering rates in 
dark matter detectors. We then introduce in section [3] four representative extension of the Stan- 
dard Model that can lead to enhanced neutrino interactions at low energy. In particular, we will 
discuss neutrino magnetic moments, light new gauge bosons with various coupling structures, and 
sterile neutrinos. The phenomenology of these models is then analyzed in section |4] with regard to 
neutrino-electron scattering and in section [5] with regard to neutrino-nucleus scattering. In sec- 
tion |6j we discuss several mechanisms that can lead to diurnal, semi-annual and annual modulation 
of the neutrino count rate at the Earth: The varying Earth-Sun distance, oscillation with very 
long (~ 1 AU) oscillation length. Earth matter effects, sterile neutrino absorption, and direction- 
dependent detection efficiencies. Finally, in section [7j we discuss existing experimental limits on 
light gauge bosons, and we show how they constrain the models introduced earlier. We also use 
Borexino and GEMMA data to set new limits on models with light force carriers. We summarize 
and conclude in section [H 



^ In a slight abuse of nomenclature, we will call the extra neutrino flavors "sterile" even though they can interact 
with ordinary matter through A' exchange. 

^ The notation keVee refers to the "electron equivalent energy in keV" , which is defined as the reconstructed recoil 
energy under the assumption that it is carried by an electron. For nuclear recoils, only part of the recoil energy 
is visible in the detector — an effect which has to be corrected for by dividing the visible energy by a quenching 
factor — so that the energy threshold for nuclear recoils is higher than that for electron recoils. When referring to 
a nuclear recoil energy, we will use the notation "keVnr" . 



3 



2. NEUTRINO INTERACTIONS IN DARK MATTER DETECTORS IN THE 

STANDARD MODEL 



Solar neutrinos may scatter elastically with nuclei or with electrons in a target and produce 
low energy recoil events. Within the Standard Model, the interaction rates are very low, beyond 
the reach of current and near future dark matter detectors, but well within the reach of the larger 
dedicated solar neutrino experiments Borexino [29j and SNO |30l I31j . Solar neutrino interactions 
will also become relevant to dark matter experiments once these experiments reach a sensitivity to 
dark matter-nucleon scattering cross sections of order 10~^^ cm^ fl]. 

The differential cross section for neutrino-electron scattering in the Standard Model is easily 
found to be 



dEr 2t:EI 



4.si{2El + El - Er{2Ey + rue)) - 2sl{Erme - 2EI) + Et 



dEr 2-kEI 



Asl{2Et + E'r - Er{2E^ + me)) + 24(^,me - 2Et) + E^ (2) 



for electron neutrinos (fg) and muon/tau neutrinos (f^, Vr), respectively. In these expressions, Gp 
is the Fermi constant, Sw is the sine of the Weinberg angle, Ei, and Er are the neutrino energy 
and the recoil energy transferred to the target electron, respectively, and nie is the electron mass. 
Note that scattering of electron neutrinos on electrons receives contributions from s-channel W 
exchange and from t-channel Z exchange, whereas for muon and tau neutrinos, only Z exchange 
is possible. The differential cross section for neutrino-nucleus scattering is 

daSM{v^.,rN ^ U^^rN) _ GlmNF'^{Er) l 2-p2 I ^7^0^2/2 IN 

dE ~ 2ttEP ^^^y^^vV^w - 1) - ErTUNSj 

+ 4Z\El + si{2El + El - Er{2E,, + m^v)) + 4(^rm^ - 2EI)))\ , (3) 

for a nucleus of mass m-Tv, nucleon number A and charge Z. The function F{Er) is the nuclear 
form factor, for which we will use the form F{Er) = 3e~'^ * /^[sin(Kr) — Kr cos(Kr)]/(«;r)^, with 
s = 1 fm, r = V-R^ - 5s2, R = 1.2A^/^ fm, k = y/2m^r (and ~ -k?) [32j. Note that for heavy 
nuclei low-energy neutrino-nucleus scattering is enhanced compared to scattering on light nuclei 
because the cross section contains terms proportional to A^, Z"^ and AZ. These terms reflect the 
fact that neutrinos scatter coherently off all the nucleons. At higher recoil energies, the form factor 
becomes relevant and compensates part or all of this enhancement. 

The maximum recoil energy of a target particle of mass rriT for fixed Ey is given by 

and, conversely, the minimum neutrino energy required for transferring a given recoil energy is 

ET' = \{^r^ ^El + 2ErmT) . (5) 

For niT ^ Er, the case we will mostly be concerned with in this paper, we can make the approx- 
imation E'™™ ~ y^mTEr/2. The event rate at a detector is obtained by folding da /dEr with the 
solar neutrino flux d^/dEu: 

dR d^ da ,^ 



dErp J ^min dEy dEf- 
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where Nt is the number of target particles in the experiment. 

In order for a nucleus of mass miy to receive a recoil energy Er ~ few keV (above the threshold 
of dark matter detectors), the neutrino must have an energy > ^ErTxif^l^ ~ few MeV. Only 
the high energy tail of the solar neutrino spectrum exceeds this threshold, so the rate of nuclear 
recoil events is low. On the other hand, the requirement for an electron recoil above threshold is 
only Ey > few x 10 keV, allowing for most of the solar neutrino flux to contribute. As a result, the 
spectrum of solar neutrino-induced electron recoil events in a dark matter detector with a threshold 
of few to 10 keV is dominated by the low-energy -pp neutrinos [Ey < 420 keV), whereas nuclear 
recoils are mostly induced by neutrinos {E^ ^ 15 MeV). 

The Standard Model neutrino-electron and neutrino-nucleus scattering rates are plotted in 
figure [T] and compared to the event rates observed in several dark matter detectors and neutrino 
experiments. Note that the neutrino-nucleus scattering rate at any given recoil energy depends 
on the target material. In figure [Tja, we have chosen germanium as an example, see [T] for plots of 
scattering rates on different nuclear targets. Note also that we neglect the fact that the electrons 
and nuclei in the target material are in bound states. This is justified because their binding 
energies are in most cases much smaller than the O(keV) recoil energies we are interested in.^ 
The only exception are the very inner electrons in the case of a heavy target material, which can 
have binding energies of order 10 keV. We have checked that including the effects of electron wave 
functions in our calculation leads only to a correction to the predicted neutrino-electron scattering 
rate of at most few-20% and introduces small spectral features at those energies where additional 
electron shells become kinematically accessible. By neglecting these small corrections, our results 
for neutrino-electron scattering become material-independent and can thus be directly applied to 
any dark matter direct detection experiment. 

In this paper we will consider a variety of new physics scenarios which add a new contribution 
to the differential cross section in the integrand of equation Obviously we will need to take 
care and respect all limits from existing measurements of the solar neutrino flux, as well as other 
neutrino experiments. It is thus useful to review the most important constraints: 

Neutrino-electron interactions: 

For neutrino-electron scattering (figure [l^) , the strongest limits on any anomalous contribution to 
the cross section comes from Borexino [29j and SNO ^30ll31j. which have energy thresholds of ~ 250 
keV and 6 MeV, respectively. This leaves little room for new physics at Er ^ 200 keV. At lower 
energies the electron recoil event rates observed in the dark matter experiments XENON- 100 |34] . 
DAMA [35] and CoGeNT [25j lie 3 to 5 orders of magnitude above the Standard Model solar 
neutrino background, so that a substantial enhancement of the neutrino-electron scattering rate 
cannot be excluded there. The most constraining measurement is by XENON-100, but we should 
stress that the precise energy threshold and rate for electron recoils below Er ^ few x 10 keV 
in XENON-100 are not precisely known since electron recoils are not relevant to XENON-lOO's 
dark matter search except as a background. The energy calibration for electron recoils depend 
on the light yield Ly, which gives the number of detected scintillation photons as a function of 
the electron recoil energy, and which was only measured at higher energies (122 keV) and has 
to be extrapolated down to lower energies. In figure [l^ we have used XENON-lOO's working 
assumption that the light yield at low energies is the same as in the calibration measurement 
at Er = 122 keV, Ly = 2.2 PE/keVee, which leads to a detection threshold for electron recoils 



^ This is difTerent for the scattering of dark matter with a mass > 1 GeV on electrons, where observable recoil 
energies can occur only when the electron enters the scattering process with a large initial momentum of at least 
a few MeV. In this case, a dark matter detector would be probing the high-momentum tail of the bound state 
electron wave functions, and, consequently, the fact that electrons are bound cannot be neglected [33) . 
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Figure 1: Expected event rates in dark matter detectors from (a) solar neutrino-electron scattering and (b) 
solar neutrino-nucleus scattering in germanium. In (a), we use units of events per keVee per year per Na 
electrons (where Na is the Avogadro number) to be able to compare rates in different materials. Thick 
black lines correspond to the total event rate, while thin lines break the rate up into contributions from 
different neutrino production processes. We also show the observed electron recoil spectra in XENON- 
100 [31] (see text for details) and Borexino , from the low-threshold analysis of CDMS data |28| , and the 
event spectra from CoGeNT [3S] and DAMA [5S|- Since CoGeNT and DAMA cannot distinguish nuclear 
recoils from electron recoils, we interpret their data as electron recoil in the left panel and as nuclear recoils 
in the right panels. 



of about 2 keVee. Measurements indicate that the light yield might actually be larger at lower 
recoil energies (down to 30 keVee) [36], and if this trend continues to even lower recoil energies, 
the energy threshold in XENON- 100 might be even lower (and the background rate per keVee 
somewhat higher) than what is shown in figure [T^. However, in many scintillators the light 
yield peaks at Er ~ 10 keVee and drops steeply below [SHI, so that the exact sensitivity of 
the XENON- 100 detector to low-energy electron recoil events remains somewhat uncertain. In 
figure [T| as well as figures [2] and [3j we indicate this uncertainty by a dashed red lines below 
Er = 50 keVee. Besides Borexino and XENON-100, also the GEMMA experiment [37j has placed 
limits on neutrino-electron scattering at low recoil energies. GEMMA limits are not directly 
comparable to to the limits shown in figure [l] because GEMMA used not solar neutrinos, but 
reactor anti- neutrinos, and thus the neutrino spectrum was different. We will comment more on 
GEMMA in section |4] when discussing neutrino magnetic moments, and also in section [7] 

Neutrino-nucleus interactions: 

For neutrino-nucleus scattering (figure ^p), we compare the Standard Model prediction to the 
observed event rates in CoGeNT |i25j and CDMS [28j. Here we discuss only elastic neutrino- 
nucleus scattering because it has been shown in [23j that the new physics contributions to the 
cross sections for inelastic processes like neutrino-induced deuteron breakup or nuclear excitations 
are about eight orders of magnitude smaller than the elastic scattering cross section in the class 
of models we are interested in this paper. The reason is that in our models, deuteron dissociation 
can only be mediated by an isoscalar vector current, whereas in the Standard Model, the isovector 
axial vector gives by far the dominant contribution. Thus, considering the 4% uncertainty of 
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the measurement [3J and the 16% uncertainty in the prediction based on the standard solar 
model [3H], no deviations from the Standard Model are expected in SNO. Note, however, that a 
dedicated Borexino search for gamma ray lines from the process ^^C + — ?■ + ^^C* — >• ffe + ^^C + 7 
could be sensitive to the types of new physics discussed here [23] ■ Also, at much higher energies 
> 100 MeV, there should be constraints from accelerator neutrino experiments such as MiniBooNE 
or MINOS. We see from figure [Ijsthat in germanium solar neutrinos can only yield neutrino-nucleus 
scattering events with Er < 10 keV because of the sharp drop-off of their spectrum at high energies. 
In materials containing lighter target nuclei, the upper end up the recoil spectrum can be higher, 
up to 20 keV for Nal and 30 keV for CaW04 [1]. Also, atmospheric neutrinos can induce scattering 
events with higher recoil energies [T] , but because of their much lower flux, we ignore these events 
in this paper. 

In the following sections, we will discuss how new contributions to neutrino-electron scattering 
and new contributions to neutrino-nucleus scattering can modify the event rates plotted in figure [T| 

3. FOUR MODELS 

We will now turn to new physics scenarios in which neutrino interactions with electrons and/or 
nuclei are enhanced to give event rates that are interesting for present and future dark matter 
direct detection experiments. As we have seen in the previous section, the precise measurements 
of the solar neutrino-electron scattering rate above few hundred keV in Borexino \29\ \39\ HO] and 
SNO \30\ [3T] set tight constraints on anomalous neutrino-electron interactions. Neutrino-nucleus 
interactions, on the other hand, should be constrained at energies > 100 MeV by accelerator 
neutrino experiments. Thus, any new physics model that could potentially contribute significantly 
to the event rate in direct detection experiments should give a signal only at low recoil energies 
around a few keV, but die away at higher recoil energies. 

The simplest possibility which fits the bill is a new interaction between neutrinos and electrons 
mediated by a very light or massless particle. In this case the matrix element for neutrino-electron 
or neutrino-nucleus scattering via the light particle contains a factor = (2mTEr)~^ , where q 
is the 4- momentum exchange, is the electron recoil energy, and rriT is the mass of the target 
electron or nucleus. This will lead to a differential cross section which is a falling function of E^. 
The new physics contribution can thus dominate the Standard Model rates at low energies relevant 
for direct detection, but not at high energies where neutrino detectors are sensitive. 

3.1. Neutrino magnetic moments 

Perhaps the phenomenologically simplest type of new physics leading to enhanced solar 
neutrino-electron scattering at low energies through the exchange of a light particle is a neutrino 
magnetic moment interaction of the form 

C^^ D /i. i^a'^^d^A^u , (7) 

where iJ,,^ is the neutrino magnetic moment and o""^ is as usual defined in terms of the Dirac 
matrices as cr"^ = IIt'^iT'^]- ^^is case the light particle which mediates new interaction is the 
photon (Aa) itself. 

While the Standard Model prediction for the loop- induced neutrino magnetic moment, fii, = 
3.2 X 10~^^//s X (mi^/eV) [411 112] (with the Bohr magneton hb = V4vra/2me, a being the fine 
structure constant) is beyond the reach of current and near future experiments, some extensions 
of the Standard Model predict sizeable fiu [HI H3tii7] . potentially close to the current 90% CL 



7 



upper limit < 0.32 x 10~^^ fiB from solar and reactor neutrino experiments |48] . in particu- 
lar GEMMA [37j.^ The differential neutrino-electron scattering rate through magnetic moment 
interactions is given by 

duaii^e — ue) 9 / 1 1 \ 

- M-Er--Er], (8) 



and the corresponding expression for neutrino-nucleus scattering is 



dE^. y Ef El, 

Here, Z is the nuclear charge, and F(Er) is the nuclear form factor (see discussion below equa- 
tion ([3])). Of course, ordinary scattering through W and Z exchange, with the cross section from 
equations ([T])-(|3]) is also present. The dependence of equation ([s]) on the neutrino energy E,y and 
the recoil energy Er arises from the interplay of the photon propagator and the derivative in the 
magnetic moment interaction vertex, equation ([T]). 



3.2. Gauged B- L 



As we have seen in section 3.1, a magnetic moment contribution to the neutrino-electron and 
neutrino-nucleus scattering cross section falls proportional to E~^ at low recoil energy. We will 
now turn our attention to scattering processes for which the recoil energy spectrum falls even more 
steeply (cx E~'^), and hence a larger enhancement of the neutrino scattering rate at low energies is 
possible without violating the Borexino constraint. 

This can be achieved if there is a new neutrino-electron or neutrino-quark interaction mediated 
by a light particle whose couplings do not contain derivatives. Let us in particular consider a 
model with gauged B — L (baryon number minus lepton number) symmetry, with the corresponding 
U{1)b~l gauge boson A' having a mass M^/ <C 1 GeV: 

Cb-l 3 -gB-L e7°^'ae + ^Qb-l qi^A'^q - Qb-l i^7"^> + ■ ■ ■ • (10) 

Here, qb-l is the U{\)b-l coupling constant and g, e and v are quark, charged lepton and neutrino 
fields, respectively. We will call A' a "dark photon" here and in the following.^ Note that we neglect 
the possibility of kinetic mixing between the dark photon and the photon here. We will discuss 
models with kinetic mixing (but with couplings to B — L) in great detail below, and we will argue 
in section [7] that, in many phenomenologically relevant processes, a coupling to i? — L is equivalent 
to kinetic mixing. 

The cross section for A'-mediated elastic scattering of a neutrino off an electron or nucleus 
depends mildly on the chiral structure of the A' couplings. Here, for concreteness we will assume 
the A' to have pure vector couplings of the form ■ijj'y^^'ipe but spectra would look similar for other 



Published astrophysical constraints can be about an order of magnitude stronger |48H50| . but it is difficult to assess 
the systematic uncertainties associated with these limits and to assign a confidence level to them. 
In the literature, the term "dark photon" is often reserved for U{1)' gauge bosons coupling to the Standard Model 
only through kinetic mixing with the photon, but we will use it in a more general context. In fact, a gauge 
boson coupled to Standard Model particles only through kinetic mixing with the photon, would not have tree level 
couplings to Standard Model neutrinos at all since they are electrically neutral. 
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chiral structures. With this assumption we obtain for the differential cross section for neutrino- 
electron scattering as a function of the recoil energy 

da A'(l^e ^ l^e) g% r o o on / ^ 

\lE. = 4.pg(Ml"2E.„..)^ [^^^ + ^' - - - -a . (11) 

The corresponding expression for neutrino-nucleus scattering is straightforwardly obtained by re- 
placing nie with the nuclear mass, and by including a coherence factor (where A is the nuclear 



mass number) and the nuclear form factor F^{Er).^ Note that in equation (11) we do not neglect 
the neutrino mass mu and we distinguish between the neutrino energy Ei, and its momentum p^, 
because later, in section [4j we will consider also scattering of heavy sterile neutrinos. 

As we will expose in section [Tj the parameter space for light U{1)b~l gauge bosons is strongly 
constrained by "fifth force" searches |52H55j which require a U{1)b-l gauge boson with couplings 
relevant to a dark matter detector to be heavier than about 100 eV. In addition, bounds on 
anomalous energy losses in stars and supernovae |56ft58j . as well as GEMMA limits on anomalous 
contributions to the scattering of reactor antineutrinos on electrons [37j severely limit the allowed 
range of coupling constants for 1 meV < Ma' ^ 100 MeV. However, we will see that some 
interesting regions of parameter space are not yet fully. Moreover, in slightly non-minimal models, 
many constraints can be easily avoided [201 EH] • 



3.3. Sterile neutrinos and a dark photon coupled through kinetic mixing 

One simple possibility to avoid many constraints on light new gauge bosons while still maintain- 
ing large neutrino-electron and neutrino-nucleus scattering rates is to consider scenarios in which 
the couplings of the new gauge boson to electrons and other Standard Model particles are much 
smaller than its coupling to at least some neutrino flavors. This is possible in models with new 
sterile neutrinos, which are singlets under the Standard Model gauge group, but charged under a 
new U{iy gauge group. Standard Model particles, on the other hand, could be coupled to the A' 
gauge boson (which we again call the "dark photon") only through a small kinetic mixing e with 
the photon. The relevant terms in the Lagrangian of the model below the electroweak scale are 



- {yhYmy^vi^ - {vsYmy^Vs - {T^LYrn^i^Us , (12) 

where the first line contains the gauge kinetic terms of the A' boson and the photon, the kinetic 
mixing term between A' and the photon with small mixing parameter e, and the kinetic term of 
the sterile neutrino I's, including its gauge coupling to the A' with coupling constant g' . (We have 
omitted the kinetic terms of the other fermions and gauge bosons in the model since they are 
not relevant to our discussion and are anyway unchanged compared to the Standard Model.) We 
could equivalently have considered the interactions of A' above the electroweak scale, in which case 
the kinetic mixing would be between the dark photon and the hypercharge boson rather than the 
photon. This would introduce not only mixing between A' and A, but also mixing between A' and 
the Z, which is, however, negligible for the low energy processes we are interested in. The second 



line of equation (12) contains the usual 3x3 Majorana mass matrix for active neutrinos ul, a 



Majorana mass matrix for the sterile neutrinos, as well as a Majorana-type mixing term. These 



^ For recoil energies below 10 keV, the inclusion of the form factor merely changes the count rate in a dark matter 
detector by 1-10% 
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mass matrices may be obtained from a seesaw mechanism if we introduce singlet right-handed 
neutrinos ur: 

1- 



- LY^HuR - UsYsH'uR - - {vrYMrvr + h.c. , (13) 

where all fields are understood as vectors in flavor space: L contains the three Standard Model 
lepton doublets, Ug contains Ug sterile neutrinos coupled to A', and ur contains S + n^ heavy right- 
handed neutrinos. We use the notation H for the Standard Model Higgs boson, H' for a new Higgs 
boson which is charged under and breaks it when it acquires a vacuum expectation value 

(vev), and we define H = e^^-fT^, where e'^^ is the totally antisymmetric tensor in two dimensions 
and a, b are SU{2)l indices.^ The Standard Model Yukawa coupling Y,^, the Yukawa coupling of 
the sterile neutrinos Yg, and the right-handed mass matrix Mr are understood to be matrices of 
size 3 X (3 + ris), Ug x (3 + n^), and (3 + n^) x (3 + n^), respectively. According to the usual seesaw 
formula, the effective (3 + ng) x (3 + Ug) Majorana mass matrix of the light neutrinos is given by 

~ \Yg {H') J ^^^R \Yg {H') J ■ 

We have not explicitly written down the kinetic and potential terms for H and H' , but we assume 
that they are such that H' acquires a small vev that gives the A' a mass consistent with the 
constraints from section [7j and the sterile neutrinos a mass that is sufficiently small to allow 
coherent mixing between active and sterile flavors at typical solar neutrino energies. Since the vev 
of H' , denoted by {H'), is by assumption much smaller than the Standard Model Higgs vev {H), 
the mostly sterile neutrino mass eigenstate will typically be even lighter than the active ones unless 
there is a large hierarchy in the right-handed Majorana mass matrix Mr. 



In models of the form (12), the sterile neutrino-electron scattering cross section is given by 
equation (11), with the replacement qb-l ^ ^/^eg' . Here, ee denotes the A' coupling to electrons, 
and g' denotes the U{1)' gauge coupling constant. We will see that in such models, daA'ive — >■ 
ve)/dEr and daA'i'^N — )• vN)/dEr can be significantly larger than in models with only the three 
active neutrinos. 

A small admixture of sterile neutrinos to the solar neutrino flux can be produced by oscillation 
before the neutrinos reach the Earth. In a two-flavor approximation with only one active neutrino 
flavor Va and one sterile neutrino flavor Ug, and with the corresponding mass eigenstates V2 and 2^4, 
the vacuum oscillation probability is given by the usual expression 

P{ya ^ ^s) = Sin' 202, sin^ (^^) , (15) 

where L is the distance traveled by the neutrinos, 624 is the effective active-sterile neutrino mixing 
angle in vacuum, and Am^g = itiI — is the splitting between the squared mass of the mostly 
sterile mass eigenstate (711-4) and the mostly active mass eigenstate (1712) in vacuum. The frac- 
tion of sterile flavors in the solar neutrino flux can be as large as 20-30% without violating SNO 
constraints on the rate of neutral current neutrino-nucleon interactions j60j . If the mass squared 
difference between active and sterile neutrinos is in the range Amli > 10~^ eV^ accessible to ter- 
restrial neutrino oscillation experiments, some of the active-sterile mixing angles are more strongly 
constrained, to the level of sin^ 26 ~ few % [8l-[10j. 

One potential problem with the model discussed here is that the sterile neutrinos (and also the 
hidden Higgs boson H' [6IJ) will acquire a tiny effective electromagnetic charge of order e [62] . 



We neglect the "Higgs portal" coupling (H^ H)(H'^ H') here. 
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There are tight hmits on such "minicharged" particles, the strongest of which come from bounds 
on anomalous energy losses in stars due to the decay of plasmons into sterile neutrinos [571 163j.^ 
Like bounds on dark photons, these constraints are also somewhat model dependent as we will 
discuss in section [7j For us, it is important to note that minicharged particles with sufficiently 
large coupling cannot leave a stellar environment and may therefore still be allowed. Also, models 
which raise the sterile neutrino mass to ~ 10 keV to several 100 keV (depending on e) completely 
evade these limits. 

Note that there are other constraints on sterile neutrinos with masses > 10 KeV, the most 
severe of which come from observations of the x-ray flux from galaxy clusters and of the cosmic 
microwave background, see |641l65j. These observables are modified if there is a large astrophysical 
population of sterile neutrinos decaying radiatively into light neutrinos. Other constraints could be 
derived from requiring the sterile neutrinos to not overdose the Universe. However, astrophysical 
and cosmological bounds can be avoided if sterile neutrinos are not produced in significant numbers 
in the early universe, for instance because of a low reheating temperature (<C 100 MeV) [66], or 
because they are chameleon-like (i.e. their effective mass could depend on the surrounding matter 
density [201 [591 [67] ) . Alternatively, the sterile neutrinos could have a fast invisible decay mode, for 
instance into light neutrinos and Majorons |64j . 

The phenomenology of > 10 keV sterile neutrinos in a direct detection experiment is essentially 
the same as that of the light sterile neutrinos discussed above, but their production in the Sun can 
be modified, especially due to kinematic suppression and because they can no longer interfere with 
the light active neutrinos. We discuss the construction of models with heavier sterile neutrinos in 



more detail in appendix ffl and their interesting phenomenology in section 4.2 



3.4. Baryonic sterile neutrinos and gauged baryon number 

While the models discussed so far are most easily detected in neutrino-electron scattering, it is 
phenomenologically interesting to consider also scenarios which predict neutrino-nucleus scattering 
to be dominant, in particular since nuclear recoils are the type of signal that most dark matter 
detectors are particularly sensitive to. A model of this type has been proposed in [23j: It extends 
the Standard Model by a gauged U{1)b (baryon number) symmetry and introduces one or several 
"baryonic" sterile neutrinos which are charged under U{1)b- Since constraints on new light gauge 
bosons coupling to quarks are much weaker than constraints on particles coupling to leptons (see 
section [T] below) , the U{1)b gauge coupling can be fairly large 



The cross section for sterile neutrino-nucleus scattering in this model is given by equation ( 11 ), 
with gB-L replaced by the U{1)b gauge coupling gB, and with the usual coherence factor and 
the nuclear form factor F'^{Er) included. 

In addition to enhancing neutrino scattering cross sections, a C/(l)s gauge boson is also a source 
of new Mikheyev-Smirnov-Wolfenstein (MSW) type neutrino matter effects |68|. 169]. In particular, 
A'-mediated coherent forward scattering of sterile neutrinos on nucleons creates a potential of the 
form 

VA' = ^{Np + Nn) (16) 



Plasmon oscillations and plasmon decay into sterile neutrinos is also an alternative production mechanism for 
sterile neutrinos in the Sun. However, the A' spectrum and hence also the sterile neutrino spectrum in this case 
are determined by thermal effects and are steeply falling functions of energy above 1 keV |:56 • We have checked 
that therefore the flux at energies accessible to dark matter experiments is negligibly small. 
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where Np and are the proton and neutron number densities in the Sun, respectively, which 
are a function of the distance from the center [38J. Note that this form of the MSW potential 



is specific to the U{1)b model. The U{1)' model discussed in section 3.3 for instance, does not 
lead to non-standard matter effects at all because the couplings of a kinetically mixed U{V)' gauge 
boson are proportional to electric charge, which is zero for ordinary matter, whereas a U{\)b-l 
model leads to a potential proportional to the number density of neutrons. (In the U{1)b-l model 
from section [3]2j this matter potential is however inconsequential for oscillation physics since it is 
flavor-diagonal among the active neutrinos.) 

The phenomenological implications of the new MSW potential are twofold: On the one hand, 
in the large potential limit 

Va' ^ max |Amj^|/2£' , j , k = 1 . . . number of neutrino mass eigenstates , (17) 

the Hamiltonian terms mixing active and sterile neutrinos are negligible compared to the large 
potential term for the sterile neutrinos. (Here, E is the neutrino energy, which is > 10 MeV for 
the neutrinos of interest to dark matter detectors.) Thus, in this case, the effective active-sterile 
mixing angles in matter are very small, and sterile neutrinos are effectively decoupled in matter. 
Their production in the Sun is then negligible, and they are only produced via vacuum oscillations 
outside the Sun. 

On the other hand, even when active-sterile oscillations are suppressed, the existence of the 
Va' term can still modify oscillations among the three active neutrino flavors. For each parameter 
point considered, one needs to check that the agreement of the theory with data from solar, 
atmospheric, reactor, and accelerator neutrino experiments is not spoiled. For solar neutrinos, 
this requires checking that the Ue — ?• I'e survival probability between ~ 200 keV-20 MeV is not 
modified substantially compared to the standard three-flavor case. Since in some models. Earth 
matter effects can lead to different survival probabilities during daytime and nighttime, one should 
also make sure that no such day-night asymmetry is present. For terrestrial neutrinos, one has 
to check that the z^e — ^ i^e and — s- survival probabilities, as well as the v fi Ve and 

— )• Vt- oscillation probabilities are not modified substantially. 

Finally, the existence of a non-standard matter potential could imply the existence of new MSW 
resonances which may lead to strong conversion of active neutrinos into sterile states. These MSW 
resonances can usually be avoided by choosing an appropriate sign for the active-sterile mass 
differences. 



4. ENHANCED NEUTRINO ELECTRON SCATTERING FROM NEW PHYSICS 

Let us now investigate the phenomenology of the models introduced in section [3] in more detail. 
We begin by studying neutrino-electron scattering rates in dark matter detectors. Most of these 
experiments make an effort to distinguish nuclear and electron recoils, focusing on the former as 
dark matter candidate events and rejecting the latter as backgrounds. Interestingly, two exceptions 
to this are DAMA [35j and CoGeNT [24j, both of which have observed a possible signal. We will 
begin with scenarios in which the scattering neutrinos (either active or sterile) are light (< 1 eV), 
and later consider also heavier sterile neutrinos. 



4.1. Scattering of Light Neutrinos 



Curve A in figure [2] shows the neutrino-electron scattering rate expected for neutrinos with a 
magnetic moment (section 3.1 ) of 0.32 x lO^^'^/zs, saturating the 90% C.L. limit from the GEMMA 
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Figure 2: Expected event spectra in a dark matter detector from new physics in the scattering of solar 
neutrinos on electrons. The different colored curves correspond to (A) a model where the neutrino has a 
magnetic dipole moment of ^i, = 0.32 x 10~^'^^s and (B, C, D) models where the scattering is enhanced by 
the exchange of a new light gauge boson A' with couplings ge to electrons and to neutrinos. The latter 
case is for instance realized in the model from section 3.3 where Standard Model particles couple to the A' 
through its kinetic mixing with the photon, but there is also a sterile neutrino Vs directly charged under U (1)'. 
To keep the discussion general, we assume the Ve transition probability to be energy-independent, and 

we have absorbed the corresponding flux suppression into a redefinition of g^. The black curve shows the 
Standard Model rate from figure [TJ and the red curves and data points show the observed electron recoil 
rates in XENON-100 [M] (see section [2] for details), Borexino CoGeNT [15], and DAMA [35]. (Note 
that CoGeNT and DAMA cannot distinguish nuclear recoils from electron recoils, so their data can be 
interpreted as either.) 



experiment [37]. We see that a significant enhancement of the event rate, by more than one order of 
magnitude at Er ~ few keV, is possible. While this is still outside the reach of existing experiments, 
near future detectors like LUX, XENON- IT, X-MASS or PANDA-X may be able to enter this 
territory. Turning this statement around, future dark matter detectors may he able to improve the 
hounds on the magnetic dipole moment of the neutrino considerably. 

Curves B, C and D in figure [2] are typical event spectra from A'-mediated neutrino-electron 
scattering in a dark matter detector. We see that, as expected, the electron recoil energy spectrum 
is proportional to the squared propagator of the light gauge boson, {q^ — M\,)~'^ where = 
—2Ej.me. It is thus a steeply falling function of for E^ > M^,/2me and flattens out for 
lower Er- This can be easily discerned by comparing curves B, C and D, which where computed 
assuming different values for Ma'- All three of these curves satisfy (but almost saturate) the 
Borexino limit, and all of them may be within the reach of LUX, XENON- IT, X-MASS, PANDA- 
X or even XENON-100, provided the detector response to electron recoils can be sufficiently well 
understood, and the electron recoil background from Standard Model processes can be sufficiently 
reduced. Even now, XENON-100 disfavors scenarios in which all of the events seen in CoGeNT 
or DAMA are explained by the scattering of solar neutrinos on electrons. It is of course still 
conceivable that only a fraction of these rates signifies the scattering of sterile neutrinos on electrons, 
and the rest is due to instrumental backgrounds. In fact, a recent preliminary investigation by 
the CoGeNT collaboration [70] suggests that such instrumental backgrounds exist. The annual 
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modulation amplitudes observed by DAMA and CoGeNT are roughly at the level of the XENON- 
100 background, and we may hope to explain these signals in models that predict a strongly 
modulating signal, see section (6) 

Curves B, C and D in figure [2] were computed with the U{1)b-l model from section 3.2 and 



the U{iy model with kinetic mixing and C/(l)'-charged sterile neutrinos from section 3.3 in mind. 
The former model is more strongly constrained since it predicts enhanced scattering rates of even 
the active neutrinos at low energies. Such signals are constrained by low-energy neutrino-electron 
scattering experiments at nuclear reactors, in particular GEMMA [37], and they may also be in 
conflict with constraints on any anomalous energy loss in the Sun and in other stars, (see section [T]). 
Astrophysical constraints can often be avoided in "chameleon" models, where the dark photon mass 
depends on the background matter density [2Ql EHl [67] . 

The U{iy model with kinetic mixing allows for much more model-building freedom than the 
U{1) B-L model as long as it is ensured that the oscillation lengths for transitions of active neutrinos 
into the more strongly interacting sterile neutrinos, 

AttE , , ( E \f eV' 
^ ~ 2.48 km X — — — 



iosc = -^=^2.48kmx — — ^ , (^ = 1,2,3), (18) 
Ami- \MeV/VAm^ ' 



are much longer than the distance of the GEMMA detector from the reactor core (13.9 m |3 
Here, we have assumed that this is the case, but that at the same time 

Lose < 1 AU (19) 
in the neutrino energy range relevant to dark matter detectors, 10 keV ^ Ey < If) MeV, so that the 



L-independent oscillatory term in equation (15) averages to 1/2 and the fraction of solar neutrinos 



converted into sterile states, P{i'e — ^ ^s), is energy-independent. Condition (19) is fulfilled for 
^Tn\^ ^ 2.5 X 10"^*^ eV^ {i = 1, 2, 3). To avoid astrophysical constraints on kinetically mixed U (1)' 
gauge bosons and on sterile neutrinos more easily, one can (but does not have to, see section [7]) 
consider sterile neutrinos heavier than ~ 10 keV or chameleon models |201 159 [ [67], 

Note that in the U{1)' model with kinetic mixing (or in any model where the signals shown in 
figure [2] originate from the scattering of sterile neutrinos), the product of couplings gi,ge needed 
to obtain the displayed curves needs to be larger in order to compensate for the mixing angle- 
suppressed flux. This is, however, easily achieved by increasing the coupling of dark photons to 
sterile neutrinos while leaving the coupling to electrons small, or even reducing it compared to 
the U{\)b~l case. In the legend of figure [2] it is understood that in models where the enhanced 
event rate is due to sterile neutrinos, the fraction of solar neutrinos converted into sterile states, 
P(fe — >■ z^s), is absorbed into a redefinition of g1. 



4.2. Heavy Sterile Neutrinos 



In the scenarios shown in figure [2j the elect ron recoil spectra in the few-10 keV region were 
falling at most as E~'^, so that a scattering rate larger than the one shown in curve B is excluded by 
Borexino, and possibly also by XENON-100 constraints. However, we note that the recoil spectrum 
has several sharp edges at higher energies of several hundred keV or above. These features are 
induced by the kinematic cutoffs of the individual processes contributing to the solar neutrino flux 
(see figure [1] for a break-down of the recoil spectrum into these contributions). We will now show 
how these kinematic edges can be shifted down to the recoil energies of interest to dark matter 
detectors. This can happen in models with heavy sterile neutrinos whose mass is close to one of 
the kinematic edges in the solar neutrino flux. Here, we will in particular focus on the sharpest of 
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these edges, coming from neutrinos produced in the reaction 

^Be + e" ^ ■^Li + Vi , (20) 

with an energy of 862 keV. 

For sterile neutrinos with a mass of that order, the velocity difference between the heavy and 
light states produced in the Sun is so large that the neutrino flux can be viewed as a completely 
incoherent mixture of heavy and light states. Hence, dynamic transitions between the active and 
sterile flavor eigenstates during propagation are absent. (There can still be coherence among 
the light states and among the heavy states, provided the mass splittings within each of the two 
sectors are sufficiently small.) We will assume that the electron recoil signal is dominated by sterile 



neutrino scattering, as is the case for instance in the 11(1)' model from section 3.3 We will also 
assume that the sterile flavor admixture to the light mass eigenstates is negligible, whereas there 
needs to be some admixture of the heavy mass eigenstate to the active flavor eigenstates in order 
to produce 1/4 in the Sun. See Appendix |A] for a discussion of models in which this is naturally 
achieved. 

Examples for the expected electron recoil rates from scattering of heavy sterile neutrinos with a 
mass not far below the ^Be line are shown in figure [sj Even though production of the heavy mass 
eigenstate in the reaction (20) is suppressed by the small mixing matrix element |C/e4p and by 
kinematic factors which are small for 2^4 masses close to the Q value of the production reaction, a 
considerable enhancement of the neutrino-electron scattering rate in dark matter experiments is 
still possible. Note that in the plot we have not taken into account kinematic modifications to solar 



neutrino production processes other than (20), many of which have 3-body final states. Instead, 
we have simply assumed the 2^4 produced in these reactions to have the same energy spectrum 
(suppressed only by the small leptonic mixing matrix element |C/e4p) as the corresponding light 
neutrinos. This introduces inaccuracies close to the upper kinematic thresholds of the pp, ^^N, 
^^O, ^^F, ^B and hep neutrino spectra, and to remind the reader of these inaccuracies we show the 
corresponding contributions to neutrino-electron scattering as dashed lines in figure [3j 

The plot shows that a model with heavy sterile neutrinos can explain an excess in low- 
threshold experiments like CoGeNT, while avoiding constraints from higher-energy detectors such 
as XENON-100 and Borexino. A certain amount of fine-tuning is, however, required in the partic- 
ular case of CoGeNT since the CoGeNT and XENON-100 thresholds are relatively close to each 
other. Such tuning of a mass splitting against an unrelated quantity (the ^Be line energy in this 
case) is reminiscent of tunings of the mass splitting against the dark matter kinetic energy in mod- 
els of inelastic dark matter [Til [72] , which were also introduced to reduce tensions among direct 
detection experiments. 

The allowed masses and mixing angles of sterile neutrinos are constrained by a large number 



of terrestrial and astrophysical experiments, as explained at the end of section 3.3 but as also 
discussed there, these constraints can be avoided in slightly non-minimal scenarios. 



5. ENHANCED NEUTRINO NUCLEUS SCATTERING FROM NEW PHYSICS 

While neutrino-electron scattering in a dark matter detector, as discussed in the previous 
section, is a very interesting discovery channel for new physics in the neutrino sector, it is not 
the process that most of these detectors are designed to look for. Let us therefore now turn our 
attention to neutrino-nucleus scattering, focusing in particular on scenarios in which the scattering 
rate at low energies is enhanced, thus possibly mimicking a dark matter signal. 

As for neutrino-electron scattering, the simplest way of achieving such enhancement is by 
introducing a neutrino magnetic moment. The expected neutrino-nucleus scattering rate for solar 
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Figure 3: Expected event spectra in a dark matter detector from A'-enhanced scattering of heavy sterile 
neutrinos on electrons (thick colored lines, see figure legend for the parameters used). We have assumed the 
A' mass to be almost negligible, and we have chosen the cross section such that the CoGeNT excess can be 
explained. Black lines show the count rate in the Standard Model, and red curves show the observed event 
rates in XENON-100 [51] (see section [2] for details), Borexino [H], CoGeNT ^ and DAMA |35|. We have 
accounted for the kinematic suppression of heavy neutrino production for the ^Be neutrinos (solid colored 
lines), but not of the pp, ^^N, ^^O, ^''F, ®B and hep neutrinos (dashed colored line), which are produced as 
parts of 3-body final states. 



neutrinos with a magnetic moment at the current upper limit is shown in figure [4j curve A, for four 
different target materials: germanium (used for instance in CoGeNT and CDMS), CaW04 (used 
in CRESST), Nal(Tl) (used for instance in DAMA), and xenon (used for instance in XENON-100, 
LUX, X-MASS, ZEPLIN, PANDA-X). As we can see, the effect is very small, and certainly not 
detectable by dark matter experiments in the foreseeable future. 

The situation is different for neutrino-nucleus scattering through the exchange of a new light 
gauge boson ("dark photon") A' with mass M^/ (see for instance the models from sections 3.2-3.4). 
In this case, the A' couplings can still be sufficiently large to allow for substantial enhancement of 
the scattering rate. Moreover, when scattering on a heavy nucleus, a low energy neutrino cannot 
resolve the nuclear substructure, and hence the scattering happens coherently on all nucleons. This 
leads to an increase in the cross section proportional to the nuclear mass number A. 

On the other hand, since nuclei are much heavier than electrons, an O(keV) nuclear recoil 
energy (above the detection threshold in a dark matter detector) requires neutrino energies of 
0(1-10 MeV), as opposed to the 0(10 keV) required for the a detectable electron recoil. This 
means that, while all solar neutrino flux components can contribute to i'-e~ scattering, only the 

and hep neutrinos — the components with the highest energy cutoff and the smallest flux — can 
affect the nuclear recoil signal. This leads to a reduction of 3 or 4 orders of magnitude when 
compared to the ^Be and pp fluxes, respectively. 

Moreover, the typical 4- momentum exchange = —2ErmN in neutrino-nucleus scattering is 
much larger than in a neutrino-electron scattering process with the same recoil energy Ej.- There- 
fore in order to obtain substantial count rates, much higher couplings are needed in comparison to 
the i^-e~ scenarios considered before. Also, the larger means that the transition between flat 
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Figure 4: Expected event spectra in dark matter detectors from solar neutrino-nucleus scattering in (a) 
germanium, (b) CaW04, (c) Nal(Tl), and (d) xenon. Note that for Nal(Tl), we use units of keVee rather 
than kcVnr for Er because due to the different quenching factors for Na (0.3) and I (0.09) [73], the nuclear 
recoil energy cannot be uniquely reconstructed. Colored curves correspond to (A) a scenario with a neutrino 
magnetic moment ^i, — 0.32 x 10^^*^/15, (B) a model with active neutrino-nucleus scattering through a light 
A' boson (for instance the U{1)b-l model from section 3.2 1, and (C), (D) a model in which 2% of the solar 
neutrino flux oscillate into a Standard Model singlet i^^, which couples to atomic nuclei for instance via a 
light U{1)b gauge boson (section 3.4) [53]. The relevant model parameters, in particular the mass of the A' 
and its coupling to nucleons {gp — gn = gn) are listed in the legend. Where applicable, we have assumed 
active-sterile mixing with sin^ 6*24 — 0.02 and Am|2 = 10^^" eV^. The black curves show the Standard 
Model rate, and the red curves and data points show the observed spectra of nuclear recoil candidates 
in CoGeNT [55], in the low-threshold data set from CDMS 28,, in CRESST i26j, and in DAMA L35j. 
The approximate Xenon-100 exclusion region is obtained by converting Xenon-lOO's observed rate of signal 
candidates (3 events/100.9 days/48 kg between 8.4 keVnr and 44.6 keVnr [751) into the units of our plots. 
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and decreasing da/dEj. happens at much larger M^/. 

In figure |4] (curves B-D), we plot the ^'-mediated neutrino-nucleus scattering rate for three 
different scenarios. The first one, curve B, involves only the three active neutrinos, assuming that 
they couple universally both to electrons and nucleons. This scenario could, for instance, be realized 
in the U{1)b-l model from section 3.2 To obtain a sizeable count rate in nuclear recoils, and not 
violate the Borexino constraint on neutrino-electron scattering, we need an A' mass heavy enough 
to suppress i^~e~ scattering, but not too heavy so that low energy neutrino-nucleus scattering is 
still enhanced. Curve B in figure [4] is a possible realization of this scenario which avoids all bounds 
to date (see discussion in section [7] and figure [s]). We see that, although the nuclear recoil rate in 
this model is higher than the Standard Model rates, it is still at least 2-3 orders of magnitude lower 
than the sensitivity of CoGeNT, CDMS, CRESST or DAMA, thus making it difficult to probe in 
present and near future experiments. 

Another possibility is to introduce sterile neutrinos which couple only to quarks, but not to 
leptons, as is the case for instance in a model with gauged baryon number and with sterile neu- 
trinos charged under it (see |23j and section 3.4). Since the active neutrinos would not feel the 
new interaction, matter effects between the active and sterile sectors would arise. Hence, oscilla- 
tion physics should be taken into account when computing the sterile neutrino flux at the Earth 
because, in principle, it could play an important role in such a scenario. To illustrate the impact 
of oscillations, we plot curve C in figure |4| for which we have assumed Am^g = 10"^*^ eV^, and 
we have chosen the couplings guQN such that the predicted event rates are comparable to the ones 
for the U{1)b~l model, curve B.^ Comparing curves B and C, we see that matter effects in this 
particular model change the rates moderately, but in an energy-dependent way. We have checked 
that oscillations of active neutrinos in solar and terrestrial scenarios are unaffected for our choice 
of parameters. 

Note that curve C could also be realized in the model from section 3.3, with a gauge group U {!)' 
under which only the sterile neutrinos are charged, and which couples to the Standard Model 
through kinetic mixing with the photon. An advantage of the U{1)b scenario is that Borexino 
limits, as well as bounds from fixed target experiments, atomic physics, and the anomalous magnetic 
moment of the muon and the electron (see section [T]) can be more easily avoided, thus opening up 
a generous window in the parameter space for mediator masses from MeV to GeV and reasonably 
large couplings. The lack of strong limits makes it even possible to have rates high enough to explain 
the current signals in CoGeNT, CRESST and DAMA. However, since only the most energ etic ^B 
neutrinos can lead to nuclear recoils above threshold in these detectors, and since the neutrino 
spectrum is steeply decreasing with energy, the event spectra cannot be fitted very well, as can be 
seen in curve D in figure El For this curve, we have again assumed Am^g ~ 10"^" eV^. 



6. ANNUAL AND DIURNAL MODULATION 

A generic prediction of almost all dark matter models is an annual modulation of the dark mat- 
ter interaction rate observed at the Earth. This modulation is caused by the relative velocity of the 
Earth with respect to the Milky Way's dark matter halo, which is larger in northern hemisphere 
summer than in winter. In fact, the DAMA collaboration has reported a statistically significant 
annual modulation in the observed event rate [271 ES] (see, however, [75l - f79] for a discussion of 
systematic effects that could cause this modulation). Recently, the CoGeNT collaboration has also 
claimed a modulating signal [25], but its statistical significance is still lower and an interpreta- 



® Our choice of the mass squared difTerence allows to have annual modulation compatible with DAMA or CoGeNT, 
see section 16.2 
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tion in terms of dark matter is problematic |80H82j (see, however, |83| which come to different 
conclusions). Both DAMA and CoGeNT observe the maximum count rate during the northern 
hemisphere summer or spring months (early June for DAMA and mid- April for CoGeNT). Other 
experiments have not yet published searches for annual modulation, but once an experiment has 
collected a sufficient number of candidate events, such a search would be the logical next step. 
Another potential "smoking gun" signature of dark matter is diurnal modulation, caused by the 
changes in relative velocity between the detector and the dark matter halo during the day. This 
modulation is predicted to be much smaller than the annual one, but it is nevertheless being 
searched for [82l[85], though with negative results so far. A significant improvement in the sensi- 
tivity to daily modulation is expected in dark matter detectors with directional sensitivity [86j. 

Because of the DAMA and CoGeNT signals, but also to scrutinize the robustness of the dark 
matter interpretation of any annual or diurnal modulation signal in the future, it is important to 
investigate alternative sources of temporal modulation in a dark matter detector. In this section we 
investigate modulation signals that could arise from neutrino physics beyond the Standard Model, 
and we discuss if and how these signals can be distinguished from a dark matter signal. 



6.1. The Earth Sun distance 



Any signal whose source is an isotropic particle flux from the Sun is expected to show annual 
modulation due to the fact that the flux decreases with the square of the Earth-Sun distance, and 
that the Earth's orbit is slightly elliptical. At its closest point, the perihelion, around January 3rd, 
the Earth-Sun distance is 0.983 AU, versus 1.017 AU at the aphelion around July 4th. Hence, 
a solar neutrino signal will be modulated with an amplitude of about 3%, with the maximum 
count rate expected in early January. The expected phase would thus be opposite to the one 
observed in DAMA, but it is important to note that, depending on the dark matter mass and 
velocity distribution, dark matter can also lead to modulating signals peaking in winter |87] . Thus, 
a neutrino signal modulating predominantly because of the ellipticity of the Earth's orbit cannot 
explain the DAMA signal, but could be confused with other dark matter signals. 

In the following, we will see that there may be additional sources of modulation that can 
potentially be large enough to overcompensate the Earth-Sun distance effect and thus reverse the 
sign of the modulation, bringing it in line with the DAMA observation. 



6.2. Annual modulation from neutrino oscillations in vacuum 



A simple way of overcompensating for the Earth-Sun distance effect is to assume that the 
anomalous signal we are interested in is exclusively due to new light sterile neutrinos (see e.g. 



the U{iy model from section 3.3 or the U{1)b model from section 3.4), and that these sterile 



neutrinos mix with the active neutrinos in such a way that the oscillation probability is increasing 
with increasing Earth-Sun distance. Then, oscillations would modify the simple scaling of 
the sterile neutrino flux with distance L, and in particular the sign of the modulation can be 
reversed [23] . This is illustrated in figure [5| where we plot the expected relative annual modulation 
amplitude, (-Rjun — Rdbc) / {Rjun + Rdbc) for sterile neutrino-electron scattering as a function of 
the electron recoil energy Er and the mass squared difference Am? between the mostly active 
and mostly sterile mass eigenstates. (Here, R stands for the count rate differential in Er, and 
the subscripts "Jun" and "Dec" indicate the count rate at the northern hemisphere summer and 
winter solstice, respectively.) Note that, in producing figure [5| we have worked in a simple two- 
flavor oscillation framework, a situation which could be realized in a 4-neutrino model, where the 
mostly sterile mass eigenstate 1^4 mixes predominantly with only one of the active mass eigenstates, 
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Figure 5: Relative annual modulation fraction (i?jun — -RDoc)/(^jun + i?Dcc) of the rate of A'-mediated sterile 
neutrino-electron scattering as a function of the recoil energy E,, and the mass squared difference between 
active and sterile neutrinos. For simplicity, we use a two-flavor vacuum oscillation framework here. One can 
clearly distinguish five different regimes, in which the rate is dominated by pp neutrinos, ''Be neutrinos, pep 
neutrinos, CNO neutrinos, and *B neutrinos, respectively. 



say z^2- We have moreover assumed that the conversion of active neutrinos into sterile neutrinos 
proceeds through vacuum oscillations. This is realized for instance in the U{iy model with kinetic 
mixing from section 3.3 in which no non-standard matter effects arise due to the zero net electric 
charge of the Sun. In models where the two- flavor /vacuum approximation is not justified, the 
phenomenology could be much richer and much more complicated. 

We see from figure [5] that a sizeable annual modulation fraction can be achieved without undue 
fine-tuning of Am? if the oscillation length Lose = ^T^Ey/ /S.w? at the relevant neutrino energies 
between few x 10 keV and ~ 15 MeV is smaller but still comparable to the Earth-Sun distance. 
Such a situation would be reminiscent of the now excluded "just-so" solution to the solar neutrino 
problem, where it was also hypothesized that the oscillation length of solar neutrinos could be 
comparable to 1 AU. It is easy to distinguish four different energy regimes in figure [5j in which the 
scattering rate is dominated by pp neutrinos, ^Be neutrinos, pep neutrinos, CNO neutrinos, and 

neutrinos, respectively. In each of these regimes, the modulation fraction is determined mostly 
by the oscillation length at the peak energy of the corresponding neutrino flux, as long as the recoil 
energy is large enough for this peak energy to be kinematically accessible. 

Note that for Am^ values larger than the ones shown in figure [s} the impact of oscillations will 
fade away once the oscillation length at the relevant neutrino energies becomes smaller than the 
diameter of the neutrino production region in the Sun. For pp neutrinos, which are produced at 
radii r < 0.2i?Q in the Sun and whose flux peaks at around 300 keV, this happens for Am^ > 
3 X 10^^ eV^, whereas for ®B neutrinos, which are produced at r < O.I-Rq and whose flux peaks 
at around 6 MeV, wash-out effects become relevant only for Am? > 10^'' eV^ [88j. 

Note also that for Am^ ~ 10~^° eV^, the Standard Model MSW potential (or the ^d'-induced 
matter potential in the case of the U{1)b models from section 3.4) will suppress the mixing of 
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electron neutrinos with sterile neutrinos at the center of the Sun. As the neutrinos propagate 
out, they can pass through an MSW resonance (if Am^ has the appropriate sign) and acquire an 
admixture of the mostly sterile mass eigenstate 1/4 because the resonance transition is non-adibatic, 
for such small Am^. (This can be seen from equation (23) below.) 



6.3. Diurnal and annual modulation from Earth matter effects 



Another mechanism by which solar neutrino signals can modulate with time is Mikheyev- 
Smirnov-Wolfenstein (MSW) type matter effects in the Earth. It is well known (see for instance 
reference [89] and references therein) that even in the standard three-flavor oscillation framework 
matter-enhanced v^,Vt z^e oscillations of solar neutrinos inside the Earth can lead to a slightly 
enhanced flux during the night, when solar neutrinos have to traverse the Earth before reaching 
a detector. This leads to diurnal modulation of the detection rate, and, since nights are longer 
in winter than in summer, it also leads to annual modulation. In the standard framework, the 
day-night asymmetry is predicted to be very small, on the few per cent level, but we will argue 
here that it can be sizeable in the new physics sector. 



Consider, for instance, a scenario based on the model from section 3.4 (a U{1)b gauge boson), 
but with two sterile neutrinos, weakly mixed with the active ones. We assume that one of the 
sterile neutrino flavors, say Vsi-, is charged under the new gauge group, whereas the other, Vs2-, is 
not. Thus, only Vgi can be observed in a detector. To simplify the discussion, we also assume both 
of the mostly sterile mass eigenstates to be heavy enough for them to be never produced in coherent 
superposition with the mostly active mass eigenstates, so that oscillations among the active flavors 
are fully decoupled from oscillations among the sterile flavors. We do, however, assume the mass 
splitting Am^ between the sterile neutrinos to be sufficiently small for oscillations among them to 
occur. We can then discuss these oscillations in a simple two-flavor framework. 

Using the well-known formalism of neutrino oscillations in matter, it is straightforward to show 
that the probability for a solar neutrino to arrive in a terrestrial detector in the Usi flavor eigenstate 
is given by (see appendix) [89] 
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Here, and are elements of the effective leptonic mixing matrix in matter at the core of 
the Sun, 9 is the vacuum mixing angle between i/gi and 1^32, E is the neutrino energy, L® is the 
distance the neutrinos travel inside the Earth, and V®, is the j4'-mediated MSW matter potential, 
equation (17), in the Earth, which by assumption affects only the flavor eigenstate Vgi^ but not 



Vs2- For simplicity, we take V^, to be constant throughout the Earth. We have also introduced the 
abbreviation 
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Note that equation (21) is valid only in the case where flavor transitions in the Sun are fully 
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adiabatic, which requires [90j 
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where V^, denotes the derivative of the non-constant matter potential inside the Sun with respect 
to distance from the center. In the following, we assume this adiabaticity condition to be fulfilled. 
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Note also that, for sufficiently heavy sterile neutrinos, V^, <^ \/^m\-\/2E, hence the active-sterile 
mixing angles in matter and vacuum are almost identical: U®^ ~ Uei, U®^ ~ C/eS- 

For oscillation lengths much shorter than the diameter of the Earth, a detector with lim- 
ited statistics and limited energy resolution will not be able to resolve the individual oscilla- 
tion peaks, but will only see their average effect. In this case we can make the replacement 
sin^(a;Am^L®/4£') — )• 1/2 in equation (21). We see that, if the resonance condition 



2EV% ~ Am^ cos 2(9 (24) 

is fulfilled, the oscillation amplitude inside the Earth can be very large even if the vacuum mixing 
angle 9 is small. If |f/^P > I^^^P, the probability for detecting a Vgi is larger during the day than it 
is at night. As discussed above, this diurnal modulation also leads to an annual modulation of the 
daily average count rate. The length of the day is larger in summer than in winter, therefore, for 
> \U%\'^ < \U®^\'^), the average daily count rate is also larger (smaller) in summer than 

in winter. Note that this type of annual modulation cannot be invoked to explain the modulation 
signals observed in DAMA |27j and CoGeNT |25j since neither of these experiments has observed 
the accompanying (and stronger) diurnal modulation [851 EI] • 

It is also very important to recall that a non-standard matter potential in the Sun (V^) and in 
the Earth (V^,) for sterile neutrinos is only generated in models in which the Sun and the Earth are 
not neutral under the gauge group that couples Standard Model particles to the sterile neutrinos. 
Therefore, the modulation mechanism discussed here will not be effective in models in which the 
new gauge boson couples to the Standard Model only through kinetic mixing. 



It is useful to rewrite the resonance condition (24) as 



ivoscVj^ = 2^cos26', (25) 

with the oscillation length Lose = 47ri?/Am^. From this expression we see that for large V®,, 
i.e. relatively strong coupling between Vgi and ordinary matter the resonance condition can only 
be fulfilled if the oscillation length is very small, i.e. for very low E or relatively large Am?. If 
-^osc ^ 1 km, matter effects are important even during daytime because detectors are typically 



located > 1 km underground. Setting Lose = 1 km in equation (25), we find that this happens for 



^>0.18GeV-2. (26) 

This should be compared to the corresponding Standard Model quantity g'^/AMyy = 8.2 x 
10-6 GeV-2. 

It is amusing to note that a strong day-night asymmetry combined with the effects of the varying 
Earth-Sun distance could conspire to give a modulation phase that is very dark matter-like. The 
amplitude of annual modulation from a day-night asymmetry is at most on the order of 20% for 
a detector at mid-latitudes, and the peak of such a modulation coincides with the summer solstice 
around June 21st. The Earth-Sun distance effect is weaker and has a nearly opposite phase, with 
a minimum on July 4th. The combined effect is a modulation with a maximum total flux occurring 
on a date earlier than June 21st and closer to early June, which is the canonical peak of the dark 
matter signal. This scenario may be easily distinguished from a dark matter signal thanks to the 
strong daily modulation which it predicts (and which has not been seen in DAMA and CoGeNT). 



6.4. Zenith angle dependence of Earth matter effects 

In the previous section, we have seen that a day-night asymmetry in the sterile neutrino in- 
teraction rate in a detector can be generated in models where the oscillation length Lose for os- 
cillations among different sterile neutrino flavors is much smaller than the diameter of the Earth 
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Figure 6: The distribution of the amount of Earth matter a solar neutrino has to travel through in winter 
(October-March, black histogram) compared to summer (April-September, red histogram) in order to reach 
a detector in the Gran Sasso Laboratory. To obtain this plot, we have computed the thickness of the rock 
around the laboratory as a function of the zenith and azimuthal angles (including effects of local topography), 
and have weighted the result by the fraction of time the detector sees the Sun in any given direction. A very 
large daily and annual modulation of a sterile neutrino signal may be achieved in models with oscillation 
lengths around a few kilometers or around several thousand kilometers. 

(~ 12 700 km), but much bigger than the detector's rock overburden (typically of order one kilome- 
ter) . Phenomenologically very interesting (though extremely fine-tuned) models can be constructed 
by choosing the mass splitting between the new neutrino mass eigenstates such that Lose is close to 
either end of this range. This is because the daily averaged distance a solar neutrino must travel 
inside rock before reaching a detector modulates annually due to the different range of zenith angles 
under which the detector sees the Sun at different times of the year. Away from the equator, the 
Sun is higher in the sky at noon in summer as compared to winter. Similarly, in winter the Sun 
drops lower below the horizon during the night. In figure [6] we show the distribution of neutrino 
path lengths in matter for solar neutrinos on their way to the Gran Sasso laboratory (where, for 
instance, the DAMA, CRESST and XENON-100 experiments are located). In the plot we distin- 
guish between the summer months (red histogram) and the winter months (black histogram). We 
see that both the average path length in matter during daytime and during nighttime varies by an 
0{1) factor between summer and winter. Thus, if Lose is around one kilometer or around several 
thousand kilometers, the oscillation probability can modulate significantly during the year (and 
also during a day). It is amusing to note that in the case Lose ~ 1 km, the precise modulation 
pattern depends on the latitude of the laboratory, on the detailed topography of the landscape 
above the detector, and possibly even on the location of the detector within the laboratory. To 
compute figure [6| we have used topographical data of the Gran Sasso region [92] to determine the 
rock thickness d{Q, ^) above the Gran Sasso Laboratory (42° 27' North, 13° 34' East, 963 m above 
sea level [93j) as a function of the zenith angle and the azimuth angle For zenith angles below 
the horizon, we have approximated the Earth as a simple sphere. We have then weighted d{Q, $) 
by the fraction of time the Sun spends at any given point (G, $) in the sky (or below the horizon) 
during the time periods indicated in the legend. All computations were done in Mathematica 8. 
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6.5. Diurnal and annual modulation from neutrino absorption in the Earth 

In models that feature sterile neutrinos with sufficiently large couplings to ordinary matter, 
the sterile neutrinos' scattering cross section can be so large that their mean free path becomes 
less than the diameter of the Earth. (For constraints on such models and a discussion of the 
allowed parameter space see section [7] below.) At night, when they have to travel through a 
substantial amount of matter before reaching a detector, the sterile neutrinos would thus loose 
all their kinetic energy and become undetectable, whereas during daytime, they could reach the 
detector unimpeded. This can lead to a very strong daily modulation of the experimental event rate, 
and due to the different length of day in summer compared to winter, also to annual modulation 
peaking in summer. 

6.6. Modulation from direction-dependent quenching factors 

Finally, there is the possibility that temporal modulation of a neutrino scattering signal is 
induced by direction-dependent solid state effects in a target crystal. This source of modulation 
is especially interesting for signals which originate from a particular direction such as the Sun, as 
opposed to signals originating from dark matter which are roughly isotropic, with only a small 
direction dependence in the velocity spectrum due to the Earth's motion with respect to the dark 
matter halo of the Milky Way. (This directionality in the dark matter velocity distribution as seen 
from the Earth is usually referred to as the "WIMP wind"). 

It is well known that the response of a solid state detector to nuclear recoils can be very sensitive 
to the direction in which the recoil nucleus is traveling with respect to the crystal axes (see for 
instance j94fi99]). In particular, if the initial momentum of the recoiling nucleus is aligned with 
one of the crystal planes, it is likely to bump into its nearest neighbors, and most of its energy will 
be converted into phonons. On the other hand, if the recoiling nucleus enters the space between 
crystal planes and travels along this "channel", it will mostly scatter on electrons, so that a larger 
fraction of the recoil energy is converted into electronic excitations. In many detectors, only 
electronic excitations can be detected, so that for these "channeled" events, the ratio between the 
visible energy and the actually deposited energy (the quenching factor) is larger. The magnitude 
of blocking and channeling effects is strongly dependent on the target material, the nuclear recoil 
energy and the temperature at which the detector is operated. In typical dark matter detectors, 
at most a few per cent of nuclear recoils with a given energy and direction can be channeled, and 
typical channeling fractions are 1-3 order of magnitude smaller |94H98j . 

Detectors like CDMS or CRESST, which use superconducting phase transition thermometers 
to measure directly the total deposited energy, are not strongly affected by channeling, and only 
their background rejection efficiency might change if the fraction of energy going into electronic 
excitations (scintillation and ionization) varies. On the other hand, experiments like DAMA or Co- 
GeNT, which rely exclusively on electronic signals for their energy measurement, could be strongly 
affected by direction-dependent quenching factors. 

In these experiments, the reconstructed recoil energy spectrum, and thus also the total number 
of events above threshold, would change as a function of the recoil direction. For conventional 
signals from the scattering of heavy dark matter particles on nuclei, the rotation of the Earth 
implies that at different times during the day, the detector sees the "WIMP wind" under different 
angles, so that one would expect a small modulation in the detection rate during the sidereal day 
(not the solar day). Direction dependent detection efficiencies would, however, not contribute to 
annual modulation because the orientation of the Earth's axis relative to the WIMP wind does not 
change during the year. This can be understood from figure [7j in which the blue cones illustrate 
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Figure 7: Illustration of the relative orientation of a detector with respect to the solar neutrino flux through- 
out the year. In blue we show the cones which the detector's normal axis traces out during a day, and in 
red we show the Earth's trajectory around the Sun. As explained in the text, a detector whose detection 
efficiency depends on the direction of the incoming particles can observe both diurnal and annual modulation 
in the solar neutrino signal, and possibly also several higher harmonics of these fundamental frequencies. 

the trajectory which the detector's normal axis traces during a day, and the red ellipse depicts the 
Earth's orbit around the Sun. 

The situation is different for signals induced by particles coming from the Sun, such as neutrinos. 
In this case, direction-dependent detection efficiencies will lead to modulation synchronized with 
the solar day (not the sidereal day) and with the solar year. The exact spectrum of modulation 
frequencies depends on the specific target material, in particular on the orientation of the preferred 
directions for channeled nuclear recoils in this material. The important time scales for modulation 
in this scenario will almost certainly be daily and annually, but semi-annual or other periods are 
also possible, depending on the symmetries and orientation of the crystal in question. 

Modulation signals from direction-dependent quenching factors may be easily identified by ro- 
tating the detector and collecting more data. 

7. THE PARAMETER SPACE FOR LIGHT GAUGE BOSONS AND NEW 

CONSTRAINTS 

One of the central ingredients of the models discussed in the previous sections is a new, light 
and feebly coupled C/(l)' gauge boson A' , the dark photon. The possible existence of such a particle 
has been considered in the literature in many contexts (unrelated to neutrinos), and strong bounds 
on the A' mass and couplings have been derived. Therefore, in order to fully assess the viability 
of our models, we need to make contact with these experimental constraints. In order to keep the 
discussion brief we will in most cases simply present the existing bounds and refer the reader to the 
literature for further details on the physics behind them. We will consider the models discussed in 
sections EjHSa 

A A U{1)b-l gauge boson with vector couplings to fermions 

B A U{iy gauge boson kinetically mixed with the photon, and with no other couplings to 
Standard Model particles 
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C AU{1)b gauge boson (gauged baryon number) 

As we have seen in sections |4] and [5| model |A] can substantially enhance even the scattering rates 
of active neutrinos. In model [BJ on the other hand, the dark photon does not couple to electrically 
neutral Standard Model particles such as neutrinos, so this model can lead to large neutrino signals 
in dark matter detectors only if, in addition to the dark photon, there are sterile neutrinos that 
carry a U{iy charge. In this case, the dark photon can induce sterile neutrino-electron and sterile 
neutrino-nucleus scattering, with the latter effect being subdominant. See sections |4] and [5] for 
details on the phenomenology. Finally, model [C] can lead to an enhancement of only the neutrino- 
nucleus scattering rate, while leaving neutrino-electron scattering largely unchanged (except for 
loop-induced effects), provided that there are "baryonic" sterile neutrinos charged under U{1)b- 

In the literature, constraints on light gauge bosons are usually presented in the context of U {!)' 
bosons coupled only through kinetic mixing (model [B]), but for many of the relevant processes 
(see below), a scenario with a U{1)b~l gauge boson is equivalent to such a model. In fact, for 
any process involving only A' couplings to electrons and protons (but not neutrons or neutrinos), 
whose electric charges happen to be identical to their B — L charges, a U{1)b-l gauge boson with 



coupling qb-l and mass Ma' is equivalent to a kinetically mixed dark photon, equation (12), with 
kinetic mixing parameter e = gs-i/e x [1 + 5^_j;^/e^]~^/^ and mass M^, x (1 + q^b-lI^j'^ ■ 
the formal level, this equivalence can be demonstrated by transforming the B — L model according 
to j4 — )• A — qb-l/g X A', and then rescaling the A' field. (Here, A is the Standard Model photon.) 
For U{1)b gauge bosons, model[C| many of the existing constraints do not apply at all (see below). 

Where applicable, we will also comment on variations of models |Af[Cl in particular scenarios with 
sterile neutrinos Vs charged under 11(1)', where A' — t- v's'^s decays can modify the phenomenology. 
We will also mention "chameleon" models, i.e. models in which the A' mass changes as a function 
of the background matter density \20\ \59\ I67j . which can also change or eliminate bounds. 

The various constraints on U{1)b-l gauge bosons are shown in the top panel of figure [sj and 
the constraints on light gauge bosons kinetically mixed with the photon are shown in the bottom 
panel. Constraints on all three models considered here are also summarized in table |lj To obtain 
the top panel of figure [sj we have translated the limits on kinetic mixing from [571 EH] iiito limits on 
a B — L gauge coupling qb-l by replacing e with gB-L/e, where e is the unit charge. (In the case 
of the constraint from T decays, the transformation is e — >• qb-l/'^g because the electric charge 
of the bottom quark is twice its S — L charge.) It is important to note that there are additional 
constraints in the U{1)b-l case (model [A]) compared to model |B] (kinetic mixing), namely from 
Borexino and from long-range force searches [5W5^ (see below for details). The latter apply also 
to the J7(1)b model (model [C]). 

The bounds in the plot are color coded: Regions colored in dark blue represent the most robust 
limits from laboratory experiments. In particular, these limits are still valid if the new light gauge 
boson is subject to strong chameleon effects. Limits in lighter shades of blue are sensitive and may 
be evaded with chameleon effects (the lightest shades of blue are reserved for limits that have not 
been worked out in detail and are less robust. Limits shown in red apply to the U{1)b-l model 
(|A|), but not models in which the A' couples to Standard Model particles only via kinetic mixing 
(model [b]). Green limits apply only if there are, in addition to the A' gauge boson, sterile neutrinos 
charged under 11(1)' with the coupling constant and mixing angle indicated in the plot. 

The individual constraints, which are also summarized in table |lj are: 

1. Anomalous magnetic moment of the electron and the muon ("((7 — 2)e", "(5 — 2)^"). 
A new gauge boson coupling to leptons (Models |A] and |B| but not [C]) will contribute to the 
magnetic moment of the electron and the muon at the one loop level, see for example |100] . 
Note that a model with parameters just below the {g — 2)^ exclusion line in figure [s] could 
explain the currently observed deviation of (5^ — 2)^ from its Standard Model prediction [101] . 
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Top: Constraints on a U{1)b-l gauge boson (model |Aj section 3.2) with coupling gs-L and 
Bottom: Constraints on hght A' gauge bosons kineticaUy mixed with the photon (model |b| 
3.3 ) as a function of the A' mass and the kinetic mixing parameter e. The various bounds are briefly 
explained in the text. Limits shown in dark blue are the ones that cannot be evaded even if the mass of the 
A' depends strongly on the local matter density (so-called chameleon effects [501 EH EH ) • The red exclusion 
regions apply only to the U{1)b-l models (|A]), but not models with only kinetic mixing (model |b|. The 
green exclusion regions apply only if the model contains sterile neutrinos directly charged under J7(l)'. Most 
limits are taken from the compilation in |57j . see text for further references. To our knowledge, the Borexino 
and GEMMA limits shown here have not been discussed before. 
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2. Fixed target experiments. Electron and proton beam dump experiments have placed 
limits on dark photons produced in the target and decaying to electron-positron pairs, 
see |102ftl04] . Obviously, this technique can only constrain dark photons with M^i > 2mf>. 
The fixed target constraints shown in figure [8] apply directly only to models |A] and [Bj In[C| 
A' decays to electron-positron pairs would be loop-suppressed, so the bounds would become 
much weaker in this case. Even in model |B| the constraints are not robust if the model 
contains sterile neutrinos directly charged under U{1)' (the case of most interest to us in 
this work), since in this case, the branching ratio for A' decay into e~'"e~ pairs is greatly 
reduced. It is expected that fixed target constraints will be significantly improved by the 
APEX experiment (see [105] for first results from this experiment). 

3. T decays. Decays of T mesons to a photon plus a dark photon, with the latter decaying 
further to /i"*"//", are constrained by i?-factory experiments jlQSj. Note that for a U{\)b-l 
or U{1)b boson these constraints are modified by an 0{1) factor compared to the case of a 
kinetically mixed dark photon because the electric charge of the bottom quark is different 
from its i? or i? — L charge. In model |b| the constraint is avoided if a fast A' decay mode 
to sterile neutrinos exists. 

4. Atomic physics constraints. By comparing the energy differences between excited atomic 
states to the energy differences measured between lower-lying atomic transitions, anomalous 
corrections to the Coulomb force at atomic distance scales can be constrained |107) . These 
constraints apply only to models in which the dark photon couples to electrons, i.e. models 
[X]and[Bj but not[C} 

5. Supernova 1987A ("SN1987A"). A dark photon may be produced in the core of a su- 
pernova and contribute to its energy loss. By requiring that in Supernova 1987 A, the energy 
loss in dark photons was not larger than the known energy loss in neutrinos (10^^ erg/s), 
constraints on the A' mass and coupling can be derived [58] ■ Note that these bounds cover 
only a limited range of A' couplings: for too small couplings, the energy loss in dark photons 
is small; for too large couplings, the supernova is opaque even to dark photons, so that 
anomalous energy losses occur only in a thin outer shell. Note also that A' emission in a 
supernova is mostly due to A' radiation off protons and neutrons. Therefore, the constraints 
apply to all three models considered here. Let us also remark that supernova constraints 
can be avoided altogether in models where the dark photon feels strong chameleon effects, 
so that its effective mass inside the supernova is higher than its mass in vacuum |20|. 

6. Solar constraints ("Sun"). A dark photon can also be produced by thermal radiation in 
the Sun. Unlike for supernovae. A' radiation in stars is dominated by emission off electrons 
or, more precisely, by the conversion of plasma excitations (so-called plasmons) into dark 
photons. Requiring that the solar luminosity in dark photons is smaller than the known 
luminosity in photons places a strong bound on A' bosons in a wide mass window around 
the solar plasmon resonance at M^i ~ 100 eV [56]. As for the supernova case, this bound is 
not robust at large e, where dark photon absorption inside the Sun becomes relevant. To our 
knowledge, dark photon dynamics in this regime has not been considered in the literature, 
but since it is clear that dark photon absorption will reduce the problematic anomalous solar 
energy loss due to A' emission, part of the region labeled capture in Sun" in figure [s] may 
still be allowed [56]. Indeed, the results from |108j (derived for keV-scale scalars) suggest 
that ©(few X 10 keV) particles may still be acceptable, while lower masses can be ruled out 
by requiring anomalous heat transfer mechanisms inside the Sun to be small. We expect 
solar constraints to be most relevant to models [A] and |B| in which the dark photon couples 



28 



to electrons. In model [C| the only available A' production mechanism is radiation off nuclei, 
which we expect to be much smaller than radiation off electrons. However, to our knowledge, 
this has never been worked out in detail. 

If the A' boson couples to neutrinos, the solar constraint also depends on neutrino dynamics. 
In the U {1)b-l model (modeljAj), where the A' boson couples directly to the active neutrinos, 
low-energy neutrinos produced in plasmon decay can carry away energy even if the dark 
photon cannot because its mass is outside the plasmon resonance region. In the literature, 
the corresponding constraints are usually referred to as "minicharged particle limits" [571l63j 
and in figure [Sj they are approximately indicated by the region labeled "Sun/Old Stars, 
energy loss via . At very large Qb-L, energy loss through low-energy neutrinos is not 
effective because these neutrinos cannot leave the Sun due to their large scattering cross 
sections. However, since neutrino scattering is suppressed by g%_i, we estimate that the 
large gs-L values required for this to happen are already disfavored by other constraints, for 
instance atomic physics constraints and g — 2 limits. Moreover, in scenarios with very large 
gB-L, also the spectrum of higher energy solar neutrino from nuclear fusion may be distorted, 
in potential conflict with experiments. While a detailed study of these issues is beyond the 
scope of this work, we estimate that in the U{1)b-l model, the region labeled "A' capture 
in Sun" in figure [s] is ruled out by minicharged particle limits. In the U{1)' model with 
kinetic mixing (model [B]), on the other hand, minicharged particle constraints are relaxed 
since the A' in this case does not couple to active neutrinos, but only to hypothetical sterile 
neutrinos. If the latter are heavier than ~ 10 keV, many minicharged particle limits no 
longer apply [63]. (For e > 10^^, even larger masses, > few hundred keV-several MeV may 
be required to avoid the constraints from [63].) Even if the sterile neutrinos are lighter, they 
may avoid solar constraints more easily if e is large. The reason is that for a 17(1)' gauge 
coupling of order one, the scattering of sterile neutrinos is suppressed only by (just like A' 
absorption, and compared to g%_i in the U (1) b-l model), i.e. sterile neutrinos in this model 
can be more easily absorbed in the Sun. For instance, for Ma' = 1 keV, e = 10~^, the sterile 
neutrino mean free path is less than 10 cm. The average energy loss per scattering is of order 
10 eV, i.e. a 1 MeV neutrino will have lost all of its kinetic energy after traveling for ~ 10 km. 
We therefore estimate that the 17(1)' model with kinetic mixing, and with M^i ~ 1 keV, 
e > 10~^, is not robustly excluded, even if it contains sterile neutrinos acting as minicharged 
particles. An in-depth study of the dark photon and sterile neutrino dynamics in the large 
e regime is a possible direction for future work. 

As a final remark, we note that chameleon effects could also help circumvent the solar energy 
loss bound if the A' effective mass is higher in the solar core than in vacuum \20\ I67j . 

Cooling of stars in globular clusters. Evolved stars, on the horizontal branch in the 
Hertzsprung-Russell (temperature vs. luminosity) diagram, can be used to set limits on 
energy loss due to dark photons in the same way in which solar constraints are derived. A 
particularly promising target for the study of these stars are globular clusters since the initial 
conditions for star formation in these objects are the same throughout the cluster. Since 
horizontal branch stars have a higher temperature than the Sun the corresponding limits 
reach out to higher A' masses. Like the constraints from solar cooling, they also disappear 
at large A' coupling > 10~®, as approximately indicated in figure [sj They also disappear in 
chameleon models, whereas their applicability to models without A' couplings to electrons 
(model [C]) has not been studied in detail yet. 



8. The CAST experiment. Helioscopes such as the CERN Axion Solar Telescope 
(CAST) [109] ■ which are looking for electromagnetic signals in a dark, shielded cavity are 
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sensitive to dark photons from the Sun entering the cavity and oscihating back into visible 
photons inside. This process can be used to set strong Umits on the dark photon couphng 
in a mass region where the expected dark photon flux from the Sun would be large |56]. 
Helioscope bounds suffer from similar model-dependencies as solar and stellar constraints; 
in addition, they can be avoided if the dark photon can decays, for instance to neutrinos, 
before reaching the Earth. 

9. Light shining through walls ("LSW"). In this type of experiments, one directs an 
intense laser beam onto an opaque wall and searches for residual signals behind the wall. 
Conversion of laser photons into dark photons and back into visible photons could lead to 
such signals, and thus these experiments can be used to set limits on the dark photon mass 
and couplings [ ITO] . These limits do not apply to U{1)b bosons (model [C|). 

10. CMB constraints. A dark photon would mix with the visible photon in a frequency 
dependent way, and this effect can attenuate the black body spectrum of the CMB 
which was measured precisely by the FIRAS instrument on board the COBE satellite |112] . 
These bounds may be avoided if the mass of the dark photon is higher in the dense early 
Universe due to chameleon effects. 

11. Borexino. Due to its large fiducial volume and the large solar neutrino flux, Borexino is 
quite sensitive to modifications of the neutrino-electron scattering rate. The non-observation 
of any anomalous signals can be translated into constraints on the couplings of dark photons 
to neutrinos and electrons. In particular, we require that the neutrino-electron scattering 
rate in Borexino should not be more than 8% [29] above the Standard Model prediction |38j . 
The solid red line in the top panel of figure [8] shows the Borexino constraint for the case of 
a U{1)b-l model (mode|A]), whereas the green lines in the bottom panel of figure [s] are for 
models of type[B](a kinetically mixed C/(l)' gauge boson), but including also sterile neutrinos 
charged under 17(1)' . 

12. GEMMA. The GEMMA spectrometer at the Kalinin Nuclear Power Plant [37] searches 
for anomalous contributions to neutrino-electron scattering, with the aim of constraining 
neutrino magnetic moments. To derive limits on the couplings of a dark photon to electrons, 
we use the data shown in figure 8 of reference |113j ^^ and compare it to the A' model using 
a simple test. We require the difference in between the A' model and the Standard 
Model to be < 3.84, corresponding to a one-sided 90% C.L. upper limit. 

Since GEMMA is looking for neutrino-electron scattering, it is insensitive to the U{1)b 
model (model [C]), and also to the U{1)' model with kinetic mixing (model [B]), unless there 
are sterile neutrinos charged under C/(l)' and heavy enough to be produced in oscillations of 



active neutrinos over distance scales of order 10 m (see section 4.1). 



13. Fifth force searches. On distance scales > 10 nm (M^/ < 100 eV), new long range forces 
are tightly constrained by precision tests of the gravitational, Casimir, and van der Waals 
forces. Since these experiments probe interactions between electrically neutral bodies, they 
are sensitive to U{1)b-l or U{1)b gauge bosons (models [A| and [C]) , but not to gauge bosons 
coupled only through kinetic mixing (model [B]). 

Note that the strong hmit on the neutrino magnetic moment derived in |113| was based on a theoretical calcula- 
tion |114| that was later revised [1151 We therefore do not consider this limit here, but we can still use the 
data presented in [113] to derive our own limits on the parameter space of dark photons. 
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Results from fifth force searches are often reported as Hmits on an anomalous Yukawa- 
type contribution to the gravitational potential (see for instance section 24 in [S3] or fig- 
ure 28 in |52j) of the form T4ew = —Giymim2/r x ac exp(— M^/r). Here, Gn is the 
gravitational constant, mi, m2 are the masses of the test bodies, r is the distance be- 
tween their centers of mass, ac is a coupling constant, and Mj^/ is the mass of the new 
force mediator, or equivalently the inverse range of the new force. Constraints on ac 
can be translated into constraints on a, B — L gauge coupling constant gs-L according 
to g%_j^ = aGGN{jni/qB-L,i){'m2/qB-L,2), with qB-L,i and qB-L,2 being the - L charges 
of the test bodies. In practice, rrij/qB-Lj — 0.5 GeV is just the nucleon mass multiplied 
by the inverse neutron fraction of the target material. (For the case of a U{1)b force, qb-l 
should be replaced by qB, and mj/qsj — 1 GeV.) 

Among the most sensitive tests of the inverse square law of gravitational interactions over 
laboratory distance scales are the experiments carried out by the Eot- Wash group [SJ [55] . 
These experiments restrict the coupling strength of aU (1) b-l boson to be smaller or similar 
to the strength of gravitational interactions for M^/ < 0.01 eV (see for instance figure 1 of 
reference |54j; for the U{1)b case, the replacement qb-l — ^ gBV^N/{Z + N) ~ O.SgB should 
be made in this figure). Such a gauge boson could thus never lead to an observable signal 
in a dark matter detector. 

At somewhat larger masses or smaller distance scales, the strongest constraints come from 
tests of the Casimir force, i.e. the small attractive force that the zero point energy of the 
electromagnetic quantum field induces between conducting objects that are brought very 
close to each other (distance 10~^-10~^ m) [52l [53]. At these distances, the boundary 
conditions imposed on the electromagnetic field by the conducting test bodies lead to a 
measurable modification of the electromagnetic vacuum energy. At even short distances, 
down to 10~^ m, also tests of van der Waals forces become important |52| [53] . Constraints 
on gB-L from tests of the Casimir effect are weaker than Eot- Wash limits for M^/ < 0.01 eV, 
but they extend to larger gauge boson masses and still provide a limit gB~L ^ 10^^ at 
Ma' = 100 eV. 

We should mention that the small allowed region between the globular cluster and fixed target 
constraints can be partially ruled out in models with only a dark photon (but no sterile neutrinos) 
using the requirement that the intergalactic diffuse photon background is not modified by A' — >• 37 
decays, that the black body spectrum of the CMB is not modified, and that the effective number 
of light degrees of freedom at the epoch of Big Bang nucleosynthesis agrees with observations |117] . 
The constraints can be avoided if the dark photon decays dominantly into neutrinos (active neu- 
trinos in the U{1)b-l model and sterile neutrinos in the 17(1)' model with kinetic mixing). 

Moreover, for dark photons with O(keV) masses or above, the dark photons themselves can act 
as "lukewarm dark matter" |117j . potentially overdosing the Universe. Constraints of this type 
can be avoided by the same mechanisms as for heavy sterile neutrinos (see discussion at the end 



of section 3.3) 



It is worthwhile to recapitulate which regions of parameter space are most relevant to the phe- 
nomenology of neutrino signals in dark matter detectors. For the kinetic mixing model (model [B]), 
dark photon masses in a window from ~ 10 keV to ~ 1 MeV, are of most interest. Masses below 
1 eV can also be considered, provided that sterile neutrinos are sufficiently heavy (> 10 keV- 
several hundred keV) to avoid bounds on minicharged particles. The kinetic mixing parameter e 
required to produce sizeable signals depends on the sterile neutrino flux and could vary between 
10~^^ and 10~^. For the U{1)b-l model (model[A| there may be a small region of interest around 
M4/ ~ 50 keV, gB-L ~ 10~^. For other masses or couplings, chameleon effects can be introduced 
to evade astrophysical and fifth force limits. The U{1)b model (model O) is far less constrained 
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"By looking for hadronic A' decays, fixed target experiments could set limits even on a ?7(1)b gauge boson (model C) . 
To our knowledge, no such analysis has, however, been done. 

'studies of hadronic T decays could be used to set limits, but to the best of our knowledge, no such limits have 
been published. 



Table I: Applicability of the various constraints on new light gauge bosons to the three models considered 
here. V indicates constraints that apply to a given model (the exact strengths of these constraints may still 
differ by an 0(1) factor between models), 'X' indicates constraints that are not applicable, and '?' stands 
for constraints for which a dedicated study would be required to determine their applicability. 

than the other two and can yield interesting phenomenology in the whole mass range considered 
(1 GeV and below). 

8. CONCLUSIONS 

In this paper, we have discussed the rich phenomenology of standard and non-standard solar 
neutrino signals in dark matter direct detection experiments. In particular, we have considered 
models featuring a neutrino magnetic moment, as well as scenarios with a "dark photon" A' (a 
light, weakly coupled new gauge boson). We have shown that in these scenarios neutrino-electron 
and neutrino-nucleus scattering can be much stronger than Standard Model weak interactions at 
the low recoil energies to which dark matter detectors are sensitive, while being consistent with 
constraints from higher energy experiments such as Borexino and SNO. If we moreover assume 
that a small fraction of solar neutrinos has oscillated into new "sterile" flavors, whose couplings to 
the dark photon are much less constrained than those of Standard Model particles, the scattering 
rates can be large enough to explain at least some of the recently reported anomalous signals from 
dark matter experiments. 

We have also discussed possible sources of temporal variations in the neutrino count rate, in 
particular the annual variation of the Earth-Sun distance (possibly in conjunction with oscillation 
effects, if the active-sterile oscillation length is not too much smaller than one astronomical unit), 
neutrino absorption in the Earth, diurnal and annual modulation due to Earth matter effects, and 
temporal modulation effects in detectors whose efficiency depends on the direction of the incoming 
particles, for instance because of ion "channeling" effects. 

We conclude that dark matter detectors, with their low energy thresholds, can be very sensitive 
to new physics in the neutrino sector, in particular if the new effects are strongest at low energy 
and are thus hidden from dedicated neutrino experiments, whose recoil energy thresholds are at 
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least a few hundred keV. Signals from neutrino-electron scattering and neutrino-nucleus scattering 
can be easily confused with dark matter scattering, especially when the signal consists of only a 
few events so that detailed investigations of the event spectrum or of temporal modulation effects 
are not possible. This is an important consideration for the present generation of detectors, but 
even more so for future experiments, whose large mass and reduced background will make them 
more sensitive to both dark matter and neutrino signals. On the other hand, our study also shows 
that dark matter detectors are powerful tools to constrain or discover neutrino physics beyond the 
Standard Model. 
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Appendix A: Models with heavy sterile neutrinos 

In section [4] we have considered heavy sterile neutrino with masses of order several hundred keV, 
which led to interesting features in the recoil energy spectrum (see figure [3]) and helped us avoid 
constraints on neutrino-electron scattering from XENON-100. On the other hand, in order to get 
a large enough scattering rate at low energies, the dark photon which mediates the interactions 
between sterile neutrinos and electrons has to be fairly light, at or below ~ 10 keV. Furthermore 
the heavy sterile neutrino needs to mix with the light active neutrinos to be produced in the Sun. 
The product of couplings and mixing angles geguSinO should be of order 10"^^, i.e. with Qu ^ 1 
from perturbativity and Qe ^ 10~^ (see figure |8|, this translates into the requirement smO > 10"^ 
Achieving this within a see-saw model similar to the one introduced in section [3] (equations (12) 



and (|13[)) is not possible because in such models the mass of the sterile neutrino is of order 
{H')'^ /Mn <^ {H') , M^i when the right-handed neutrinos are very heavy. 

Apart from the neutrino masses, we also need to generate an appropriate mixing matrix. In 
particular, in order for the sharp drop in the spectra visible in figure [3] to occur, the electron recoil 
signal has to be dominated by scattering of the heavy mass eigenstate rather than by scattering 
of the light mass eigenstates i^i, f2, through their small admixture. (Remember that it is 
the the sterile flavor eigenstate Vg which is charged under the new gauge group.) This requires 
the leptonic mixing matrix element C/e4 (which gives the 1/4 admixture to the electron neutrino) to 
be much larger than the elements Ugi, Us2, Us3 (which give the lyg component of the light mass 
eigenstates). 

Building a model that satisfies all these requirements requires some engineering, but is certainly 
not impossible, as we now show. We add a pair of sterile neutrinos, one of them (say, i'sl) charged 
under C/(l)' and the other one (z^^^) neutral (anomalies can be taken care of by introducing extra 
spectator fermions). The Dirac mass of the sterile neutrinos is thus proportional to {H'), which 
we will take to be around 800 keV. The mass of the dark photon is also proportional to {H'), but 
may be smaller if the 17(1)' gauge coupling is smaller than one. For instance, let us assume a gauge 
coupling of order 0.01 to give the dark photon a mass around 10 keV. The presence of the H' vev 



33 



also allows us to write mixing terms between active and sterile neutrinos, particularly between the 
left handed sterile and the right handed active states. The full neutrino mass matrix is 
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Y* (H) 











Mr 














y; {H') 






Mil, 




Mr^44 


where we have placed all neutrinos in a 


vector ■01/ = 







(Al) 



and flavor indices are 

suppressed. If there are, in addition to the three active neutrino flavors and their right-handed 
partners, left-handed sterile neutrinos and the same number of right-handed sterile neutrinos, 
Yi, and Mr are understood to be 3 x 3 matrices, Y,^i, Yg, and M/j 14 are 3 x matrices, and Yg, 
Mr^44^ are Ug x Ug matrices. 

To obtain a model with the desired masses and mixing angles, we will assume that 1^4, Mr^h, 
and Mij^44 are negligibly small. This could, for instance, be justified by declaring sterile lepton 
number (carried by I'gL and i'sr) a global symmetry only broken by the spurion y,'. Another 
possibility is to introduce two C/(l)'-breaking Higgs fields, H[ and H'2, with only H[ carrying 



on 



sterile lepton number. Then, the term Yg* {H') in equation (Al) would be replaced by Yg* {H[ 
and the term Y* {H') would be replaced by Y* {H2). The terms Mr^u, Mr^^^, and Y*^^ (J 
the other hand, would be absent. 

If we set Yiy4, Mr^i4, and Mr^44 to zero and assume Mr ^ (H) , {H') in equation (Al), we find 



that to leading order in the small masses, the admixture of the 0{{H')) mass eigenstates to the 
left-handed active flavor eigenstates is 



|C/e4| 



Yuy;, {h) 

V2Yg Mr 



(A2) 



The admixture of the C/(l)'-charged flavor eigenstate Ugi to the light 0{{H) /Mr) mass eigen- 
states, on the other hand, vanishes to leading order as desired. 

10^^ GeV will produce 
100 GeV can 



As is clear from equation ( A2 ), a standard seesaw mechanism with Mr 
mixing angels which are too small. However an electroweak-scale seesaw with Mr 
easily yield mixings that are large enough. For example, taking Yg and Yg of order one, and 
requiring the active neutrino masses niy ~ {Yy {H))M^^{Yy (H)) to be at the 0.1 eV scale, gives a 
mixing angle of 



\U, 



e4 



10" 



100 GeV\ 
Mr J ■ 



(A3) 



In order to produce spectra that are close to those observed by DAMA and CoGeNT one can pick 
gu 1 and ge ~ 10~^, which is not ruled out by the constraints from section [t} 



Appendix B: Earth matter effects on solar neutrinos 



In this appendix we explain in detail how Earth matter effects in the sterile neutrino sector can 
lead to daily modulation (see section 6.3 ) |89] . We assume a model with two sterile flavor eigenstates 
I'gi and i'g2, with i^gi giving the strongest signal in the detector. Thus, we are interested in the 
probability for z^e — ^ i^si transitions. As in section 6.3, we make the simplifying assumption that 



the new, mostly sterile, mass eigenstates and 1^5 are so heavy that they cannot be produced 
in a coherent superposition with the light mass eigenstates 1^1, 1^2, t's- We do assume, however. 
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that coherence between and 1^5 is possible due to a small mass splitting. As the neutrinos are 
produced in the core of the Sun, and as we assume the oscillation length to be much smaller than 
the size of the Sun's core, the solar neutrino flux at the Earth can be treated as an incoherent 
mixture of z^4 and 1^5. Hence, the computation of the probability fg — )• V31 is two-step: first, we 
need to determine the amount of z^4^5 that exit the Sun; then we need to calculate the probability 
of detecting these mass eigenstates as Vgi after they have traveled a distance L® in the Earth. 
Finally, we perform an incoherent sum, that is we compute 

P{Ve -^Vsl) = Peil^e iyi)P®il^4. J^sl) + Peij^e 1^5) P® {^5 ^ l^sl), (Bl) 

where P0 and P^ are the transition probabilities in the Sun and in the Earth, respectively. 

Since only the two flavors Vsi and z^s2 are relevant to us, we can define an effective two-flavor 
rotation matrix with mixing angle 6 

U'=( ''''^] . (B2) 

\ — sin cos 6 J 

We denote the vectors of the vacuum mass basis, flavor basis and the effective mass basis in Earth 
matter by \i^a), and respectively. The relations between them are given by 

\u^) = U'jW^) and \u^) = U^*\uf), (B3) 




where U® is the effective mixing matrix in Earth matter, which has the same form as equation ( |B2[ ), 
but with 9 replaced by the effective mixing angle in Earth matter, 6^, which is obtained when the 
full Hamiltonian is diagonalized. Concretely, working in the flavor basis we have 

Am2 J \ J 

(with Am^ = ml — ml) and 

jjm HU® = diag(-Amy4S , Am^/4^) . (B5) 

It is easy to show that 

Am2cos20-2SF® ^ . Am2sin20 

cos 2(9® = — ^ ^ and sm26'® = — ^ , (B6) 

Ami Am| 

where 



Am| = uj ^m^ = Am^y sin^ 26 + (cos2 29 - 2EV®//\m'^Y . (B7) 

Now, we calculate -P©(z^4 — )• i^si) = Ki^si|e~*^^®|z^4)p. With the abbreviations cg = cos 9, se = 
sin 9, c© = COS0© and s© = sin^©, we flnd 



|(i/,i|e-^^^®|i/4^'' 




g , iAm|Lg3/4£; 




-iAm2 Lg,/4£; 

(B8) 

cos^ 9 - sin^ 29 ' sin^ ( | (B9) 



.2 a ..„2o^2i?F^ /u^Am^Lg 
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In the first line we liave used the relations (B3 ), and the third line then follows from simple algebraic 
manipulations. It is straightforward to see that to obtain i-©(j^5 — )• I'si) we just need to replace 
cos^ 9 by sin^ 9 and reverse the sign of the second term. 

To complete the calculation of P(fe — )• z^si), we need Pq^i^s ^i,^)- The probability for 
producing the matter eigenstates corresponding to v^^^ at the center of the Sun is given by \Uf^ gp, 



and since the adiabaticity condition (23) is well fulfilled (unless sin 20 is tiny, in which case we 



would not expect an observable signal from sterile neutrino scattering anyway), this is also the 
probability that a solar neutrino exits the Sun as i^4^5, i.e. 



ir/® P 

I^e4,5l ' 



(BIO) 



Plugging (BIO) and (B9) into (Bl), we obtain equation (21). 
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